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Abstract 

We deduce a new formula for the perihehon advance of a test 
particle in the Schwarzschild black hole by applying a newly developed 
non-linear transformation within the Schwarzschild space-time. By 
this transformation we are able to apply the well-known formula valid 
in the weak-field approximation near infinity also to trajectories in the 
strong-field regime near the horizon of the black hole. The resulting 
formula has the structure = ci — C2 ln(c3 — e^) with positive constants 
ci,2,3 depending on the angular momentum of the test particle. It is 
especially useful for orbits with large eccentricities e < C3 < 1 showing 
that — )• 00 as e — )• C3. 

Keyword(s): Perihelion precession, perihelion advance, Schwarzschild black 
hole 



1 Introduction 

Surprisingly little is known about the value of the perihelion advance in the 
strong-field region of the Schwarzschild space-time, especially at large eccen- 
tricities of the closed orbit. Recent interest in calculations and measurements 
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of orbital characteristics like the relativistic pericenter precession, also called 
relativistic perihelion advance, see e.g. [1], [2], [3] and [4], show the necessity 
to have good closed-form expressions for it in the most prominent example 
of a space-time, the Schwarzschild black hole. Also the calculation of perihe- 
lion precession in classical mechanics is subject to recent interest, see [5] and 
section 5 of [6] . For latest work on analogous questions in axially symmetric 
space-times see e.g. [7]. In section VII of [8], the nearly circular motion of 
a particle within /(i?)-gravity has been discussed, whereas in [9], the orbital 
motion in /(i?)-gravity with quadrupole radiation was calculated. 

It is the aim of the present paper to apply the ideas and formulas deduced 
in section 10 of [6] to the calculation of the perihelion advance of periodic 
orbits in the Schwarzschild black hole.^ Especially, we want to find out, 
how B depends on the eccentricity e of the orbits. To have a stricter posed 
question, we ask: what is the change of O if we change the eccentricity e of 
the orbit but keep the angular momentum h = r^ip of the particle^ constant? 

The core of the deduction is a newly developed non-linear transforma- 
tion within the Schwarzschild space-time. Due to its importance we have 
chosen to present this transformation in two independent versions, one by 
the set of eqs. (2.16) to (2.36), which directly applies to the geodesies in 
the Schwarzschild space-time, the other one is given in the appendix, which 
is a self-contained and more abstract deduction of the primarily unexpected 
symmetry of eq. (2.16). 

The notation is as follows: Let the orbit be the periodic but non-constant 
function r{(fi) with invariantly defined Schwarzschild radius r, then © is de- 
fined by the period 27r -|- © of the function r, i.e. r{ip) — r{(p -|- 27r -|- ©). The 
eccentricity e is defined by 

e='-^ (1.1) 
r2 + ri 

where r2 = ma.xr{ip) and ri = min r{ip). Thus, e and are invariantly 



^Of course, the exact formula is well- known: it contains elliptic integrals, but in practice 
this formula is of minor use only. And, probably more importantly: these elliptic integrals 
which can be evaluated numerically to every degree of accuracy, do not easily lead to 
the identification of the physically interesting quantities we are trying to find out in the 
present paper. 

^The dot denotes the derivative with respect to the eigentime of the particle. 
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defined quantities, even for orbits which may be very far from being eUiptic 
ones. Moreover, e does not change if r is replaced by c • r with a constant c. 
If we identify e and — e, then also the replacement of r by its inverse 1/r in 
eq. (1.1) leads to the same eccentricity. 

By continuously changing ri and r2 we get circular orbits at ri — r2, 
and interchanging ri and r2 leads to the same geometry of the orbit, hence 
the same value of G. This means, we expect to be an even function 
of e at constant angular momentum. From now on we restrict e to the 
interval < e < 1, nevertheless, this consideration is useful, as the Taylor 
development of around e = should only contain even powers of e. 

The reason why we parametrize the periodic non-circular orbits by an- 
gular momentum is the following one: At fixed value h, these orbits can be 
uniquely be parametrized by e, but can also be uniquely be parametrized by 
the perihelion ri of the orbits, and can also be uniquely parametrized by the 
total energy E of the test particle. This is not a trivial statement, as for the 
general case, given h and more than one orbit exists: e.g. one hyperbolic 
orbit and another one leading towards the horizon. 

2 Geodesies in the Schwarzschild black hole 

Wc take the Schwarzschild solution in Schwarzschild coordinates^ with mass 
parameter m > as usual: 

ds^^(l-—\ de - r^dn^ (2.1) 

V r J l-2m/r ^ ' 

where dil^ is the metric of the standard 2-sphere. We apply units such that 
fight velocity c—1 and Newton's gravitational constant G —1.'^ To calculate 
time-like geodesies we may assume without loss of generality that they are 
situated in the equatorial plane. The angular coordinate is denoted by </? and 
the proper time along the time-like geodesic (^(t), r(r), '^{t)) is denoted by 

^But see [10] for clarifying historical notes to this notion. 

*In principle, we could also apply units such that G = 1/m, and under these circum- 
stances, we have one less parameter in all the calculations, but then the departure from 
the usual well-known formulas is becoming even larger. 
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T, a dot denotes d/dr. We assume t > 0. Then we get from eq. (2.1) 

2m\ -2 2-2 



1 r - — — - rV ■ (2.2) 

We restrict to the region outside the horizon, i.e. to r > 2m. 

h = rV (2.3) 

is the conserved angular momentum of the test particle. We exclude purely 
radial motion which is characterized by /i = and choose the orientation of 
space such that /i > 0. Inserting eq. (2.3) into eq. (2.2) we get 

2m\ P 

1 . 2.4) 

r J 1 — 2m/r H 

A further conserved quantity is the energy E defined by 

E^(l-'^)i>Q. (2.5) 



r / 

Inserting eq. (2.5) into eq. (2.4) we get 

1 = -, . 2.6 

To remove f, the third of the three terms containing a r-derivative from eq. 
(2.6), we describe the path of the particle as r{(p) and get via ^ = ^ by the 
help of eq. (2.3) 

h dr 

Inserting eq. (2.7) into eq. (2.6) we get after multiplication with 1 — 2m/r 

To get complete information from the geodesic equation, we still need the 
radial part of it, we take it from eq. (9.6) of [6]: 

1 - 2m/r r3 (1 - 2m/r)2 ' ^ ' ' 
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Following [11], we introduce the nondimensionalized inverted Schwarzschild 
radius u via 

u = -— (2.10) 
m ■ r 

and the dimensionless parameter e via 

.-^>0. (2.11) 

This leads io r = h? / {m ■ u) and 

dr du 
dip m ■ dip 

Inserting eq. (2.12) into eq. (2.8) we get 



With eq. (2.11) we finally get 



ev?'\ f_ 2eu\ _,2 e (du^'^ 



(2.12) 



We multiply by 3/ {2e) and get 

1 ( du\'^ V? eu^ 3/^9 \ /^-.^x 

- — -u + = — (£;2-l). (2.15) 

2\dip) 2 3 ^ 2£ V ) ^ ' 

Derivating this equation we get 

d'^u ^ 2 



+ u^l + su\ (2.16) 



dip^ 

Eqs. (2.15)7(2.16) represent the motion of a particle u in the potential 

V{u) = -u + ^-'-^ (2.17) 
with fjL interpreted as energy. A dash denoting ^ we get 

V {u) = -1 + u - eu'^ (2.18) 

and 

V"(u) ^l-2su. (2.19) 
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To shift the turning point of this potential to the origin, we define 

v = u-^ (2.20) 
as new variable instead of u. Eqs. (2.17)/(2.20) lead to 

11 1J EV^ 

Defining further 

then eq. (2.15) now reads 

I ( dvV /I \ ev^ , , 

o +^hr-i -^ = /^2. (2.23) 



2 \d^) V4£ 

Let us note the special solution f = 0for£ = l/4 and /i2 = 0. It corresponds 
to a special semistable circular orbit. Apart from this special solution, eq. 
(2.23) possesses periodic solutions v{ip) only for the parameter range < £ < 
1/4. This we will always assume in the following. 
With the notation 

V,{v) = ^ - l) - ^ (2-24) 

we get 

Vl{v) = ^-l-sv'^ (2.25) 

and 

V^'iv) = -2ev . (2.26) 
The equation ('^) — possesses the solutions 

±^ (2.27) 

representing one maximum and one minimum. This forces us to introduce 
the new variable 

w = '^'"^ (2.28) 

instead of v. With 

.3^i^ (2.2.) 
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we get now from eq. (2.23) 

In the next step we replace ip hy a new angular coordinate 

= v^l -4£. (2.31) 

It is to be observed, that ^ = is identified with (p — 27r, so = is 
identified with ^ = 27r • ^1 - 4s. Eqs. (2.30)/(2.31) together with 

/.4 = ^ (2.32) 

lead to ^ 

1 f dw\ w , . 

This is the equation we are going to solve now. Derivating eq. (2.33) we 
arrive at the simple equation 



dijj'^ 2 ■ 
Remarkably enough, it contains no e. The function 

fH = f - y (2-34) 
has zeroes at w = and w — ±-\/3. For the derivative we get 

fH = ^ - Y (2-35) 

possessing zeroes ai w = ±1. /(—I) = —1/3 is the local minimum and 
/(I) = 1/3 is the local maximum of /. We note that also /(— 2) = 1/3. The 
constant solutions of eq. (2.33) are w{tp) = —1 with /i4 = —1/3 representing 
the stable circular orbits and w(^) = 1 with /X4 = 1/3 representing the 
unstable circular orbits. Besides these exceptions it holds: Every periodic 
solution w{ip) of eq. (2.33) is a non-constant one and is completely confined 
in the interval — 2 < w{ip) < 1. This is related to the energy parameter 114 
being confined to —1/3 < /i^ < 1/3. Let wi — maxw{ip) and W2 — minu'('0). 
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then //4 = f{wi) — f{w2). We have — 1 < < 1 and —2<W2< —1. Let 
us parametrize these periodic orbits by the parameter 64 defined by 



Wl = 264-1 with < 64 < 1 . (2.36) 

In the hmit 64 — > we get the stable circular orbits. 

In the other limit 64 — > 1 we come arbitrarily close to the unstable circular 
orbits, that means, that the perihelion advance tends to infinity in this limit. 

Using eqs. (2.34)/(2.36) we get 

Me,) = f{w,) = -le^ + 26^ - 1 . (2.37) 

Due to dfi^/de, = 464(1 — 64) > this represents a one-to-one correspondence 
between 64 and /i4. Solving now the equation ^4 = f{w2) for W2 we get with 
eq. (2.37) 

«^2 = ^-64- ^3-^1- ^-64) . (2.38) 

Let ■00 = k{e4) be the complete period of the function w{ip). We calculate it 
by solving eq. (2.33) with w{0) — W2 and w{ipo/2) — wi and get 

^0 = ^(64) = 2 / , . (2.39) 

-^"^2 ^2/^4 ■u;3/3 - w 

Though the exact solutions for such integrals can be found in the hterature, 
see e.g. [12], page 355, for the general theory and [13], [14] for its concrete 
application, they are of little use as the elliptic integrals can be evaluated 
by the Weierstrass function only, and not in the usual closed-form presen- 
tation, which would allow for a physical interpretation. The three cases 
64 e {1/2, 0, 1} will now be considered in detail. 

Let 64 = 1/2, then fj,4 — wi — and W2 — —y/S. Eq. (2.39) simplifies to 

/O fjiij 
-^/3 . A„3 /Q _ 



-V3 ^w^/2,-W 

The substitution x — y/—w / leads to 

A;(l/2) = 4- ^-7 (2.41) 
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where 







according to standard tables, the Gamma-function has r(l/4) = 3.62561, 
and we get / = 1.31103 and 

k{l/2) = -00 = 6.90164 = 27r • 1.09843. (2.43) 

To get higher accuracy, the 1.09843 in eq. (2.43) has to be replaced by a 
with 

The hmit 64 — > corresponds to — > —1, and /"(— 1) = 1, so the system 
represents a harmonic oscillator with unit frequency: 

^/;(^) = -1 + 64 -cos V', A;(0) = Vo = 27r . (2.44) 

The limit 64 — ?> 1 can similarly be solved: it corresponds to Wi — )■ 1, and 
/"(I) = —1, so wc have to replace cos by cosh in eq. (2.44) and 64 by 1 — 64: 

^(■0) = 1 — (1 — 64) • coshz/) , ■^o — >■ 00 ^ 64 — >■ 1 . (2-45) 

To quantify the diverging part we replace cosh^^ by exp('0)/2 and 1 — 64 by 
(1 - e|)/2, then eq. (2.45) reads 

w(^) = l-exp(t/')(l-e^)/4. 

w{'^) — will be reached at = lii(4/ (1 — e|)), so the general structure 
must be approximately 

k{ei) = ^0 = ci - C2 • In (1 - e^) (2.46) 

with certain positive constants ci and C2 of order 1. We will fix them by the 
condition that eq. (2.46) holds exactly true for 64 = 1/2 and 64 = according 
to eqs. (2.43) and (2.44) resp. This leads to Ci = 27r and C2 = 27r • 0.34215. 
So eq. (2.46) reads 

k{e^) = 27r • (1 - 0.34215 • In (l - e^)) . (2.47) 
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For 64 <^ 1 we get approximately 

k{e4) = 27r • (1 + 0.34215 • el) . (2.48) 

To get higher accuracy, the 0.34215 in eqs. (2.47), (2.48), (2.51) and (2.52) 
has to be replaced by (3 with a from above and 

a — 1 



In4-ln3 



Let us now return from the ■(/'-picture to the (/9-picture. Let (po be the 
complete period in the 99-picture, then we get by eq. (2.31) 

^0 = = . (2.49) 

The perihelion advance is © = (/Jq — 27r, i.e. 

e = -^k^ - 27r . (2.50) 

Inserting eq. (2.47) into eq. (2.50) we get 

This eq. (2.51) is useful for all 64 with < 64 < 1 and a strict result for 

64 = 1/2; also both the limiting behaviours 64 ^ and 64 ^ 1 represent 
strict ones. 

Inserting the approximation eq. (2.48) into eq. (2.50) we get 

This eq. (2.52) is useful for all 64 with < 64 <^ 1. Inserting the expression 
(2.11) for £, we get the final formula, expressing the perihelion advance in de- 
pendence on angular momentum h and the parameter 64, which is correlated 
to the eccentricity of the orbit: 

^^^^ 1^1 -0^34215. to 
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3 Properties of the orbits 



To be able to interpret eq. (2.53) we need a better knowledge of the orbits. 
The formula (1.1) for the eccentricity is invariant with respect to multipli- 
cation of r by a constant factor, and also (after identification of e with — e) 
invariant with respect to a replacement of r by 1/r, but not invariant with 
respect to adding a constant to r. But just this addition of the constant has 
been done in eq. (2.20), so it must be expected, that the parameter 64 will 
depend both in e and on h, and not only on e. 

Due to eq. (2.11), the condition 2m < r < 00 reads < u < 3/{2e), 
and we restrict to the interval < £ < 1/4 because only for these £- values 
non-circular periodic orbits exist. The limit £ — > is the Newtonian limit; for 
this case we know the non-circular periodic orbits to be exact elliptic ones, 
i.e. © = for all values of h and e. For e — 1/4, only one periodic orbit 
exists, it is a circular one. For e > 1/4, the angular momentum is too small 
to allow for periodic orbits, the particle goes to r — > 00 or to r < 2m after 
sufficient long time. 

To find the real properties of the perihelion advance, we have therefore to 
restart at the point just before the translation (2.20) has been apphed. The 
function V{u), eq. (2.17), has the following zeroes: V{0) — 0, for e > 3/16 
this is the only one, for £ < 3/16 the other zeroes of V are calculated via 



where < uq < 3/(4£) < U5 < 3/{2e), and in the special case e — 3/16 we 
have the double root U5 — uq — 3/(4£). 

Let Ui be the maximal value of u{(p) and U2 be the minimal one. Then 
we get ^ 

m{l±e) 

^1,2 = = /q , 2^ y^-^) 

mri^2 P — f^y^j + e j 

leading to 

U1 + U2 ^ m 

2 p - m(3 + e2) ^ > 



^The parameter p has its usual meaning, its relation to ri,2 and e can be seen e.g. from 
eq. (3.9) below. 
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and 

ui - U2 ^ me 

2 p-m(3 + e2) ■ ^ ' ^ 

Solutions with constant value u fulfil u = 1 + eu^, see eq. (2.18), i.e. the 
function dV{u)/du has the two zeroes: 

1^3,4 = ^ • (l ± VT^) (3.5) 

where < -^4 < 1/(2^) < u-i < 1/e. The second derivative of V{u)^ see 
eq. (2.19), vanishes for u = l/(2e) only. Therefore, at -u = -^4 we have a 
minimum of V{u), and at m = M3 a maximum. For the special case e = 3/16 
we have = U5 = Uq = 4 and V^(4) = 0. 
We get 

V{u,) = ^^.{l-6e-{l-Asf') (3.6) 

and have obviously always V{u4) < 0. To find the sign of V{u3) more calcu- 
lations are necessary. We get 

n«3) = Y^-(l-65 + (l-4e)3/^) (3.7) 

which is positive for < £ < 3/16 and negative for 3/16 < e < 1/4. Thus, to 
get the set of non-constant periodic bounded orbits, we have to distinguish 
three cases:^ First case: For < £ < 3/16 such orbits exist for < —/j, < 
—V{u4) with n taken from eq. (2.15). Second case: For e = 3/16 such orbits 
exist for < < 16/27. Third case: For 3/16 < £ < 1/4 such orbits exist 
for -V{u3) < -jj, < -V{u4). 

We parametrize the bounded solutions by ri and r2, where ri is the 
perihelion and r2 the aphelion. In the present calculations we restrict to the 
parameter values 2m < ri < r2 < 00, i.e. to motion completely outside the 
horizon r = 2m. 

We pose as additional restriction to the possible values of ri and r2 the 
property, that a bounded orbit with these values really exists. 

Second case: s = 3/16 implies h = 4m. All values < —fi < 16/27 are 
possible, leading to 25/27 < E'^ < 1, with E from eqs. (2.5)/(2.14) and all 

^The second case can be subsumed to the first one by allowing e < 3/16 as well as 
under the third one by allowing e > 3/16. But as the second case possesses other special 
properties it is simpler to deal with it in an extra case. 
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values of the eccentricity < e < 1 are possible. In the xi-picture we get: All 
values < W2 < 4/3 an all values 4/3 < iti < 4 are possible; and this also 
exhausts the set of all possible values. In the r-picture we get: perihelion ri 
has 4m < ri < 12m, and aphelion r2 has r2 > 12m. In the limit e — )■ we 
get the perihelion advance 6 27r(-\/2 — 1) and in the limit e ^ 1 we get 
the perihelion advance O — )> oo. 

First case: < £ < 3/16, i.e. h > 4m, and all values of the eccentricity 
< e < 1 are possible. Physically, this first case can be interpreted as 
follows: Further increasing the energy of the particle at constant angular 
momentum would have the consequence that its path goes to r — )■ oo, so 
no periodic orbit appears. This behaviour we know already from Euclidean 
geometry: if the eccentricities of a set of ellipses tend to 1, then the result is 
a parabola. 

Third case: 3/16 < e < 1/4, i.e. 2v^m < h < 4m, and all values of 
the eccentricity < e < e(e) are possible, where e{e) < 1 is the following 
expression: 

e(£) = (3.8) 

In the limit e — > 3/16 we get, as expected, e{e) — > 1. It is, however, quite 
surprising, that in the other limit £ —> 1/4 we get e{e) — > 0. The perihelion 
advance tends to infinity if the eccentricity tends to e{e). For £- values being 
only slightly below 1/4, the perihelion advance becomes quite large even for 
extremely small eccentricities. Physically, this third case can be interpreted 
as follows: Further increasing the energy of the particle at constant angular 
momentum would have the consequence that its path goes to r < 2m inside 
the horizon, so no periodic orbit appears. 

To get an easy comparison with other deductions we also define the arith- 
metic mean a = (ri + r2)/2 of ri and r2 and call it semimajor axis. The 
geometric mean of them is denoted by 6 = y^rv^r^ and we call it semimi- 
nor axis. Even if the orbit is not an exact ellipse, we use the usual formula 
(1.1) for defining the eccentricity e = (r2 — ri)/(r2 + ri) having the range 
< e < 1. Further, we use the parameter p = b'^/a which is also sometimes 
used in deahng with ellipses, it holds p = a(l — e^) = 2rir2/ (ri + r2), so p is 
just the harmonic mean of ri and r2- 
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For the exact ellipse in the Euclidean plane we have the equation 

P 



r — - 

1 + e cos if 

thus 

and r2 = -^. (3.9) 
1 + e 1 — e 

Therefore, the parametrization of the set of orbits discussed here with (ri, 

can also be done with (p, e), and the latter one eases the comparison with 

the literature. 

In the next step we prescribe the values of ri and r2 and calculate the 
energy E and the angular momentum /i > by use of eq. (2.8). At the 
extrema, dr/dip — 0, so we get as one of the conditions 

Inserting eqs. (3.9) into eq. (3.10) we arrive at 



= /i^ (p - m(3 + e^)) (3.11) 



which makes sense for p > m(3 + e^) only. This inequality will always be 
assumed to hold in the following, it is equivalent to 

n>m--^, (3.12) 

meaning that the minimally allowed value for ri depends on the eccentricity, 
it holds 

3 + , , 

2m < m ■ < 3m . (3.13) 

1 + e ^ ^ 

Analogously we get 

3 + 

r2> m > 3m . (3-14) 

1 — e 

Inequality (3.12) can also be expressed as ab'^ > m(4a^ — 6^). Prom eq. (3.11) 
we get 

h= , (3.15) 
y'p - m(3 + e^) 

and then 

E^ ^ (3.16) 
^Jp{p - m(3 + e^)) 
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Further, the parameter s calculates to 

_ _ 3m ?>m^{?> + e^) 

One easily calculates 



(3.17) 



and for p — 2m(3 + e^) one has equality in the middle of this chain of 
inequalities. 



4 Discussion 

Let us interpret eq. (2.53), i.e. 

g^^^ 71 -0.34215. to (1- el) _ 

V ^1 - 12mV/i2 

This fourth root in the denominator seems at least a httle bit dubious, so 
we compare with the literature. First we concentrate on the orbits with 
neghgible eccentricity but allow strong fields, i.e., the orbit may be close to 
the horizon. Then eq. (4.1) reduces to 

e = 27r-|^ ^ -1 I ^27r- I - 1 | . (4.2) 



Second, we apply 1/v^l — 5 = 1 + 5/4 + 55^/32 + . . . and get the weak- field 
limit by inclusion of the first two terms to 

^ ^ /3m2 45mn 
Eq. (4.40) from [14] reads in our notation 

As one can see, eqs. (4.3) and (4.4) become identical in the hmit e — > 0, 
thus confirming the correctness of our calculation at least in this order of 
approximation. 
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To have a better comparison with other known results, we continue the 
discussion with the exact circular orbits. In our notation they can be calcu- 
lated by inserting r — into eqs. (2.6) and (2.9), i.e.. 



1 + 



(4.5) 



and 



J.3 _ 2m/ry 

resp. These equations can be solved for h and E by 



(4.6) 



h = ^ and E = ' . (4.7) 

Jl — 3m/r Jl — 3m/r 

Therefore, only for r > 3m such orbits are possible, and in the limit r 3m, 
a light-like circular orbit appears: The velocity of the particle is 



\ mir 

^ ' (4.8) 



1 — 2m/r 

which tends to 1 as r — )■ 3m. For every r > 3m, a circular orbit exists, and 
with h from eq. (4.7) we calculate the parameter e — "irv? jli? to 

^-^V (4.9) 



e— — . ^ , 
r V r / 

If we insert the expression e eq. (4.9) from the circular orbits into eq. (4.2) 
we get 

e = 27r • I ^ ^ - 1 I . (4.10) 

V^l - 12m(l -3m/r)/r 

This is exactly the same as 

e = 27r ■ I ^ ^ - 1 I (4.11) 



-y^l — 6m/7 



which represents the expression for orbits close to circular ones already de- 
duced in [6], eq. (11.2). For m <^ r we develop eq. (4.11) to 

^ 67rm 67rm 
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a form which can be found in the majority of texts, e.g. as eq. (1) of [13]. 
In the weak-field approximation apphed to eq. (4.8) we get velocity ^Jm/r, 
and for the Tinit mass test particle E'km = m/(2r), E^^t — —fn/r, hence total 
energy E = 1 — m/ (2r). 

Let us finally discuss the case 64 = 1/2 in more details. Then the perihe- 
lion advance reads exactly also for strong fields, see eqs. (2.43)/(2.53): 

3^^^ / 1.09843 _\ 

V^l- 12mV/i2 J 

We want to calculate the properties of the orbits to which this case belongs, 
the eccentricity should be something like e = 1/2, but, as already said, the 
relation between e and 64 may slightly depend on h. Looking back to eq. 
(2.40) we see that fi^ = and Wi = maxw = and 'W2 — minw = 
With e = 3m^/h^ and eqs. (2.29)/(2.32) we have //g = //s = 0. With eq. 
(2.28) we then get: 

vi — maxv — and V2 — minv — — VS — 12e/2e . 

With eq. (2.20) we then get: 

ui = maxti — — and U2 = minw = — • f 1 — \/3 — 12^) . 
2e 2e ^ 

To represent a closed orbit, 1^2 > is necessary, i.e. l/6<£<l/4. Applying 
r = /i^/ {mu) leads to the perihelion 

ri = 6m (4.14) 

and aphelion 

r2 = 6m/ (1 - V3 - I2e) . (4.15) 
This gives rise to the eccentricity 

^ ~ 2/V3- 12£^ ^^'"^^^ 

thus having e — > and — > 00 as £ — > 1/4; and e — > 1 and © — > 27r • 0.4456 
as £ —> 1/6. This means: For all periodic orbits with perihelion r = 6m, eq. 
(4.13) represents an exact formula for the perihelion advance. 
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To find out, how the perihehon advance changes with eccentricity, we 
solve eq. (4.16) for e: 



1 1 



4 3 (e + l)2 
Then eq. (4.13) reads 



e = 27r ■ ■ 1.09843 • ^Jl + 1/e - ij , (4.17) 



where 1.09843 is an abbreviation for 

This shows that at least in this range, a Taylor development around e = is 
not possible. 

Next, we want to find out, how the perihelion advance changes with 
aphelion 

r2 = 6m ■ (4.18) 

leading to 

e = 27r • (^v^ • 1.09843/-\/l - 6m/r2 - 1^ . (4.19) 

The formulas deduced here allow to evaluate orbital properties, especially 
the perihelion advance for all periodic orbits in the Schwarzschild field, and 
in regions near the horizon they are much better and easier to handle than 
other methods used in the literature. The method used has the potential 
also to be apphed to alternative theories of gravitation, see e.g. [15] and [16] 
and the references cited there, to get the theoretical background for possible 
experimental tests of the theories. 



5 Appendix 

The main idea of this paper is outlined by the set of transformations between 
eqs. (2.16) and (2.36). To find out the internal structure of this idea, we 
present in this appendix another and independent deduction of the unex- 
pected non- linear transformation of eq. (2.16), i.e. 

<P"^ -.2 /r 1X 

+ ^ = i + (5-1) 
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which was the core of the calculations in the main part of this paper. It is 
a self-contained deduction, so it can be read independently from the main 
text, and shows the internal symmetry of eq. (2.16)/(5.1). 

Here we concentrate solely on this eq. (5.1) and how solutions at different 
values e are related to each other. Here, e < 1/4 is a positive real parameter 
and u = u{(p) is a non-constant periodic function of the independent variable 
(p. Let the smallest period of m be -(9 > 0. Due to the simple structure of 
the equation one can easily show that i) can also be defined as follows: If 
u{ipi) is a local minimum and u{Lp2) is the next one, then u{(pi) = u{(p2) and 
= ^2 — ^i- Thus, the perihelion shift for this orbit is 

e = '(?-27r. (5.2) 

Now we fix a further real parameter A > 1 and define the function v = v{(p) 
as follows: 

v{ip) = B - uiVXip) - A (5.3) 

where A and B depend on e and A only, and -B > 0. Clearly, f (v?) is 
also a non-constant periodic function. The function v has smallest period 
■& = i^/VX and perihelion shift & = '& — 27r. With eq. (5.2) we get 

e = e/VA-27r(l-l/VA). (5.4) 

Thus, the perihelion shift of the function v is smaller than simply the perihe- 
lion shift of the function u divided by \/A as one would have expected from 
a first glance. In this sense the transformation discussed here is a non-hnear 
one. Now we want to fix A and B such that 

— + v^l + ev^ 5.5 

with a parameter e = e{e, A). This means, that v shall solve essentially the 
same equation as u does, only the value of the parameter e may differ. To 
find the form of A, B and e we insert (5.1) and (5.3) into (5.5). This leads 
via 

(fv _ ^ cPu{^/~\(p) 

to 

B\[l + su^ - + Bu{-\/\ip) - ^ = 1 + e[B'^u^ + - 2ABu] . (5.6) 
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Here, u means u(-\/X(p). This equation (5.6) must be identically fulfilled. The 
vanishing of the terms proportional to leads to BXe — B'^e, i.e. to 

(5.7) 

Inserting (5.7) into (5.6) wc get 

EX - B\u{V\^p) + Bu{V\^p) - A = 1 + A^\e/B - 2A\eu{V\^p) . (5.8) 
The vanishing of the terms proportional to u leads to 

-BX + B ^ -2AXe , 

i.e. to 

One gets A> 0. Finally, inserting (5.9) into (5.8) we get 



2Xs B \ 2Xe J ' 

To ensure S > we restrict to the region where A^(l — 4^) < 1, i.e. where 

1<A<^^=. (5.11) 
Inserting (5.10) into (5.7) and (5.9) we get 

e^\-[l-X\l-Ae)]>Q (5.12) 



and 

2(A-1) 



1 - A2(l - As) 

resp. To find the upper limit for i. we calculate 



(5.13) 



e — e 



(i-^)-(A^-l) (5.14) 



which leads to < £ < £. 
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The limit £ — > is a singular one, as in this limit, the non-linear equation 
(5.5) becomes a hnear one. Nevertheless, it leads to the following result: All 
non-constant solutions of eq. (5.5) with e — have the period ■»? = 27r, i.e. 
perihelion shift = 0. With eq. (5.4) we get for this case 

e = 27r(yX-i) . (5.15) 

With eq. (5.12) we get for £ = 

A = (1 - 4£)-^/' . (5.16) 

Combining eqs. (5.15) and (5.16) we finally get for the perihelion shift 

This is exactly the same as eq. (4.2). 
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